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ABSTRACT 

I present the Standard Model calculation of the decay rate T(b —* sg) (g 
denotes a gluon) at next-to- leading logarithms (NLL). In order to get a 
meaningful physical result, the decay b — > sgg and certain contributions 
of b — > sff (where / are the light quark flavours u, d and s) have to be 
included as well. Numerically we get £> NLL (6 — > sg) = (5.0 ± 1.0) ■ 10~ 3 
which is more than a factor 2 larger than the leading logarithmic result 
B LL (b — ► sg) = (2.2 ± 0.8) ■ 10~ 3 . Further, I consider the impact of this 
contribution on the charmless hadronic branching ratio £y, which could 
be used to extract the CKM-ratio IT^/V^I with more accuracy. Finally, 
I have a short look at £y in scenarios where the Wilson coefficient C% is 
enhanced by new physics. 
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Abstract. 

I present the Standard Model calculation of the decay rate T(b — > sg) (g denotes 
a gluon) at next-to-leading logarithms (NLL). In order to get a meaningful physical 
result, the decay b — > s<?g and certain contributions of 6 — > s// (where / are the 
light quark flavours it, d and s) have to be included as well. Numerically we get 
i3 NLL (£> — > sg) = (5.0± 1.0) • 10 -3 which is more than a factor 2 larger than the leading 
logarithmic result B hL (b — > sg) — (2.2 ± 0.8) • 10~ 3 . Further, I consider the impact of 
this contribution on the charmless hadronic branching ratio B^, which could be used 
to extract the CKM-ratio |T4b/14&| with more accuracy. Finally, I have a short look 
at in scenarios where the Wilson coefficient Cs is enhanced by new physics. 



1. Introduction 

Theoretical studies for inclusive 5-decays have become accessible due to heavy quark 
effective theory (HQET). HQET |2|,[| states that the amplitudes of decaying 5-mesons 
can be expanded in powers of (AQ CD /mft) where the leading term is nothing but the 
decay of the underlying 6-quark. Corrections to this start at 0(AQ CD /mf) only, which 
numerically amounts to ~ 5%. 

One interesting subclass of inclusive B-decays are the charmless ones B — > Xj. 
They have been thoroughly investigated in the literature: the calculation for b —>■ q'q'q 
with q = d,s and q' = u,d,s have been available to NLL for quite some time @J§J|,0] 
and the missing piece, b —* sg, is what is presented here. Due to the sensitivity of the 
charmless branching ratio £y to the poorly known CKM ratio |V^&/Vd,|, these decays 
might be used in order to get a better determination of this standard model parameter. 

The decay b —>■ sg turned out to be an interesting issue in the discussion of the 
'missing charm puzzle'; for a long time there was a discrepancy between theory and 
experiment for the average charm multiplicity n c per 5-decay and for the inclusive 
semileptonic branching ratio B a \ [§]. This issue seems to have settled a bit. 

\ Work done in collaboration with Christoph Greub Q| . Work partially supported by Schweizerischer 
Nationalfonds. 
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The calculation is based on an effective Hamiltonian 7i e ff obtained by integrating out 
the heavy degrees of freedom of the standard model (i.e. the t-quark and the H^-boson). 
Retaining operators up to dimension six, the 5-flavour effective Hamiltonian responsible 
for b — > sg reads 

AC 8 

H cS = -F-VZV^C^O^), (1) 

where are CKM-matrix entries, Cj(//) the Wilson coefficients and Oj(/z) the relevant 
operators. As the Wilson coefficients C*3_6 are small, we only consider the operators 

O x = {s L l,T A c L ) (c L YT A b L ) , 2 = (s Ll ,c L ) (c L Yb L ) , 

° 8 = leb ("^T A b R ) G A V . (2) 

Here T A stand for the SU(3) co \ om generators. The small CKM matrix element V u b as 
well as the s-quark mass are neglected. The whole work is done in the NDR scheme, 
i. e. with anticommuting 75 and using MS subtraction. 



3. Contributions 



3.1. Virtual Corrections to 0\ and O2 

In figure [l] there is the complete set of Feynman diagrams that contribute to b — > sg. 
As these start at 0(g s a s ) only (no one loop contributions), the Wilson coefficients C\ 
and C2 are needed to leading logarithmic precision only. 



Figure 1. Complete list of two-loop Feynman diagrams for b — > sg associated with 
the operators 0\ and (9 2 . The crosses denote the possible locations where the gluon 
is emitted. 

b sb sb sbs 




The heart of the procedure is to make a Mellin-Barnes representation of the 
two-loop integrals. That is, after the Feynman parametrization and after integrating 
over the loop momenta, the remaining expression is of the form (k 2 — M 2 )~ x which can 
be written as 

1 1 



(k 2 - M 2 ) x {k 2 ) x T(A) 



h) 27i / M-M 2 /k 2 yT(-s)T(\ + s) , (3) 
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where, in the actual calculation, M 2 (resp. k 2 ) is m 2 (resp. m 2 ) times some combination 
of the Feynman parameters and 7 is a path in the s-plane parallel to the imaginary axis 
with —A < 3?(s) < 0. This representation leads naturally to an expansion in z = m 2 /ml 
which is numerically around 0.1. All contributions up to z 3 are retained. In the end, 
the result consists of powers and logarithms of z only. 

The matrix elements for the operators O1/2 contain UV-divergences which are 
removed by renormalization; there are no IR problems and also singularities due to 
m s = are absent. 

3.2. Virtual Corrections to Og 

Some comments on the contributions of Og are in order. The relevant diagrams are 
displayed in figure |2|. 



Figure 2. Complete list of Feynman diagrams associated with the operator Og. The 
analogues of a), c) and d), where the virtual gluon hits the 6-quark instead of the 
s-quark, are not shown explicitly. The real gluon can be attached to any of the crosses 
shown in a). 

b 08 s b g> 08 <g> s b 08 s b g s 3 b 08 s 



b) c) gP- d) 3 e) 

■3 «3 



Some technical details of the calculation: for this part, the s-quark mass has not 
been set to zero but was kept gulator in order to distinguish between mass- 

singularities and IR-singularities. 

As Og contributes to b — > sg at tree-level already, the Wilson coefficient Cg is 
needed at NLL precision |H|. Also due to the appearance of the gluon field strength 



in the operator, Z 3 enters upon renormalization. From the gluon self-energy there are 
logarithms ln(/^/m/) (with / = u,d,s) entering the calculation. In order to avoid the 
singularities that would show up when the respective masses are set to zero, also the 
decay b — > sff (arising from 0$ only) was included. 

After renormalization Tg(6 — > sg + sff) has the form 

T 8 (b -> sg + sff) = 22 a nm (4) 

e IR 

n,m=0,l,2 
n+m<2 

where it is anticipated that the remaining singularities will cancel against the gluon 
bremsstrahlung. 
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3. 3. Bremsstrahlung b — > sgg 

In order to get a physically meaningful (in particular an IR finite) result, the 
bremsstrahlung contributions form Oi/ 2 /8 have to be added. Again (for 0$ only, as 
the ones from 0\ii are finite anyway) m s is kept gulator and the phase space 

integrals are done in d = 4 — 2e dimensions. Only the singularities are worked out 
analytically, the phase space for the finite parts is done numerically. 

4. Result 

Adding the various pieces, the singularities cancel in fact (according to the KLN- 
theorem). The result can be cast into a form involving an effective matrix element 
D (where T NLL (b -> sg) denotes r(6 -> sg(g)) + T 8 (b -> sff)): 



r NLL (& - sg ) = ^G 2 F \v*v tb \ 2 \D\ 2 + rST" (5) 

D = C 8 °' eff + g - f C 8 °' cff + C?[ii ln(m b //i) + n] 

+ C 2 °[£ 2 ln(m 6 //i) + r 2 ] + C 8 °' eff [(4 + 8 + O ) ln(m 6 ///) + r 8 ]} (6) 

where r^ ems denotes the finite bremsstrahlung part which affects the result by ~ 5% 
only; £ 1 / 2 and (£$ + 8 + (3q) are nothing but the entries of the anomalous dimension 
matrix 7i 8 / 2 8/88 an< ^ r i/2/8 encode the loop functions. For the Wilson coefficients we 
have used the perturbative decomposition Cf = Cf + f^C}. 

The branching ratio is given by the expression Bib — > sg) = - — ^^- B^. 



with 



si 



Figure 3. Scale dependence of the function D (equation (g)) and the branching 
ratio in various approximations: the dashed line corresponds to leading order (for D 
this is simply Cg' e (/*)), the solid line is the full NLL result. For the dotted line we 
set r\ = r-i = r§ = (and, for the branching ratio, rfj" oms = in addition). The 
dot-dashed line is obtained by putting ri — only. 
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The NLL branching ratio still suffers from a big /^-dependence. As is illustrated in 
figure |3| this is due to r 2 which is also responsible for the sizable enhancement of the 
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branching ratio. This large factor is multiplied by a s (/x)C2(/i) and numerically spoils 
the cancellation of the \i dependence which is established through £1/2/8 ■ 
Numerically we get (and reproduce the LL result UTTf ) 

B NLL {b^ sg) = (5.0 ±1.0) ■ 1(T 3 , B LL {b^ sg) = (2.2 ± 0.8) • lO" 3 . (7) 

Combining our NLL result with the work of Lenz et al. , we calculate the CP- 
averaged charmless hadronic branching ratio £y and the CP-averaged total charmless 
branching ratio B(B — > no charm) (the errors are estimated by varying the scale [i in 
the range [7715/2, 2mJ) 

B 4 = (1.88^°)%, B{B -> no charm) = (2.22l°j°)%. (8) 

The latter is very sensitive to |Ki&/Vrf,|: varying this ratio in the range 0.06 — 0.13 the 
total charmless branching ratio covers the range (1.50 — 2.42)%. 



5. Remarks on n c and B s \ 

As discussed in the introduction, the theoretical predictions for both, the charm 
multiplicity n c and the semileptonic branching ratio £> s i used to be in disagreement 
with the experimental data [[F]. This discrepancy decreased a lot after the inclusion of 
the complete NLL corrections to b —>■ cuq and b — > ccq (q = d, s) [[1^] . If one allows the 
renormalization scale fi to be as low as mb/4, then there is a marginal overlap between 



theory and CLEO- and LEP data [|H| . It is therefore a matter of taste if one considers 
this problem to be solved or if one is inclined towards an enhancement of the Wilson 
coefficient through non standard model physics. The impact of such an enhancement is 
displayed in figure f|. 



Figure 4. Charmless hadronic branching ratio B</ as a function of the ratio 
/ = Cf(m w )/Cf' SM {m w )- The dotted (solid) curve includes the LL (NLL) 
approximation for b — > sg. 
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f=Ci«(m w )/C^M(m w ) 
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